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Based on the systematic investigation of the data available for A > 40, a 



Z 1 ' 3 dependence for the nuclear charge radii is shown to be superior to the 



generally accepted A 1 ' 3 law. A delicate scattering of data around R c jZ x ' 3 



is infered as owing to the isospin effect and a linear dependence of Rc/Z 1 ' 3 



on N/Z ( or (N — Z)/2 ) is found. This inference is well supported by the 

h : 

microscopic Relativistic Continuum Hartree-Bogoliubov (RCHB) calculation 

conducted for the proton magic Ca, Ni, Zr, Sn and Pb isotopes including the 

exotic nuclei close to the neutron drip line. With the linear isospin dependence 

provided by the data and RCHB theory, a new isospin dependent Z 1 ' 3 formula 

for the nuclear charge radii is proposed. 
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Nuclear radius is one of the most fundamental bulk properties of an atomic nucleus 
|1],|2|]. Among all the size quantities describing nucleus, nuclear charge radii has been inves- 
tigated by various techniques and methods experimentally |3|-|9|], including the muonic atom 
spectroscopy [Q, isotope shift of optical and K X-ray spectroscopy |§-|8| and high energy 
elastic electron scattering ||, etc.. Recently more and more nuclei far from the /3-stability 
line become accessible experimentally thanks to the development of radioactive ion beam 



facilities fl0|jTl|| . The nuclear size connected with exotic phenomena such as skin and halo 
have become a hot topic. The understanding of its property has importance not only in 
nuclear physics, but also in other fields such as astrophysics and atomic physics, etc. With 
its accuracy, the study of the nuclear charge radii is very important to understand not only 
the proton distribution inside the nucleus but also the halo and skin. Particularly if one can 
get a simple and reliable formula for nuclear charge radii, it will be very useful to extract 
the de-coupling of proton and neutron in the exotic nuclei and provide information for the 
effective nucleon-nucleon interaction widely used in all the nuclear models. Here in this let- 
ter the available experimental charge radii data for A > 40 will be examined and its global 
behavior will be studied. Instead of the widely accepted A 1 / 3 law, a new Z 1 / 3 formula with 
isospin effect will be proposed. 

Based on the consideration of the nuclear saturation property, nuclear charge radii R c 
are usually described by the A 1 / 3 law [|T]j2| 

Rc = taA 1 ' 3 , (1) 

where A is the mass number and R c = */| (r 2 ) 1 / 2 , with (r 2 ) 1 / 2 the root-mean square (rms) 
charge radius. For very light nuclei, because of their small A and large fluctuation in charge 
distribution due to the shell effect with short period, it seems that the charge distribution 
radius as a bulk property has little meaning. A detail analysis of charge radius data for 
A > 40 shows that ta is by no means a constant, but systematically decreases with A; i.e., 
ta ~ 1.31 fm for light nuclei (A ~ 40) and ta ~ 1-20 fm for very heavy nuclei (see upper 
left panel in Fig. 1). This fact implies that some physics is missing in Eq. ([J). 



A definite evidence of the violation of A 1 / 3 law is also found in the measurements of 
isotope shift in mean square charge radii JT^JT^]. In particular, 5{t 2 )a+2,a values (associated 



with an addition of two neutrons) are often found to be considerably smaller compared to 
what is expected from the A 1 / 3 law (5{t 2 )a+2,a = jA~( r2 )A)- A typical example is that the 
observed charge radii of the calcium isotopes 40 ~ 50 Ca remain almost the same (except a 
very little change induced by deformation or shell effect), though the mass number A has 
changed significantly. In contrast, there is also evidence that the observed 8{t 2 )a+2,a values 
(associated with the addition of two protons) are often greater than what is expected from 
the A 1 ' 3 law (e.g., 5(r 2 ) for 46 Ti- 44 Ca, 50 Ti- 48 Ca, etc.). 

Along the /^-stability line, the ratio Z/A gradually decreases with A, i.e., for light 
nuclei Z/A « 1/2, and for the heaviest /3-stable nucleus 23 |U, (Z/A) 1/3 « 0.7285, thus 
(l/2) 1 / 3 /(Z/A) 1//3 « 1.09, which is very close to the ta ratio 1.30/1.20 shown in upper left 
panel of Fig. 1. A naive point of view is that the charge radius of a nucleus may be more di- 
rectly related to its charge number Z, rather than its mass number A. Therefore, compared 
to the A 1 / 3 law, a Z 1 ! 3 dependence for nuclear charge radii may be more reasonable 

R c = r z Z 1 ' 3 . (2) 



as noted in Ref. ||12| . An analysis of the very limited data of charge radii then available 
showed that rz remains almost a constant, i.e., rz = 1.65(2) fm for A > 40. The Z 1 ! 3 
dependence of nuclear charge radii was also used to modify the Coulomb energy term in 



the semi-empirical nuclear mass formula [ 14]] , and it was found that the agreement between 



the calculated and experimental results was improved. Moreover, the A" 1 ' 3 law for the 
nuclear giant ( monopole, dipole and quadrupole ) resonance energy (oc 1/R) also could be 
improved, if the A' 1 / 3 dependence is replaced by a Z^ 1 ^ 3 dependence [|i~5}|. 

In the past two decades, a vast amount of new experimental information on the elec- 
tromagnetic structure of nuclear ground states of many nuclei has become available |3] [| . 
and accuracy has been improved. In particular the muon factories at Los Alamos (LAMPF) 
and at Villigen (PSI, formerly SIN) started their operation at 1974. Almost all stable nuclei 



have been measured by the muonic X-ray transition technique and the corresponding charge 
radii have been rather accurately deduced (the experimental relative error is about 10~ 3 ). 
Moreover, modern techniques for optical isotope shift measurements have made it possible 
to reach even short-lived (down to 1 s) unstable isotopes ||. Therefore, it is worthwhile to 
reexamine the fundamental property of nuclei and to investigate whether the vast amount of 
improved experimental results follow the Z 1 ! 3 dependence. The values of measured (r 2 ) 1 / 2 
for 536 nuclei with A > 40 compiled in Ref. [|3|-|9| are analyzed in Fig. 1 by using the 
A 1 / 3 and Z 1 / 3 dependence, respectively. The dependence of charge radii on the quadrupole 
deformation j3 has been taken into account for the rare-earth deformed nuclei, i.e. |J 

r A = r Ad (l + -^[3 2 ), r z = r Zi {l + £-0 2 ), (3) 

and for spherical nuclei (/3 = 0): r^ = r^d, r z = r Z d, and the values of (3 are taken from 
Ref- |,0. 

In the upper left and right panels of Fig.l, the charge radii for the most stable 159 
nuclei with A > 40 along the /3-stability line have been analyzed by using the A 1 / 3 and Z 1 ! 3 
dependence. In the middle left and right panels, the same has been done for the measured 
(r 2 ) 1 / 2 for 536 nuclei with A > 40. Two significant features can be observed: (A) On the one 
hand, the agreement between the data and the calculated results using the Z 1 ^ 3 dependence 
is much better than that using the A 1 / 3 law, i.e., while there exists a global regular decrease 
of rAd with A, r Z d nearly remains constant (r Z d = 1.631(11) fm). The relative rms deviations 
a for the Z 1 ! 3 dependence [p = 7. hi x 10~ 3 for stable nuclei and 1.00 x 10~ 2 for 536 nuclei 
) are much less than that for the A 1 / 3 law (a = 1.90 x 10~ 2 for stable nuclei and 1.63 x 10~ 2 
for 536 nuclei ). (B) On the other hand, though the rms deviation for the Z 1//3 dependence 
is significantly reduced, an isospin induced scattering of the data in the middle panels in 
Fig. 1 can be also observed compared with that in the top panels. In fact, r Z d generally 
increases with iV for most isotopic chains, e.g. for ^" 96 Zr, |- 100 Mo, ^" 104 Ru, ^ 2 " 110 Pd, 

ir ll6 Cd, 112-m Sn; 122-130 Te; 124-136 Xe) 142-148^ 144-154^ 154-160^ ^ (except for 

only a few lighter isotopic chains, e.g., 3g~ 86 Kr, 3g _88 Sr, and a small anomalous decrease of 



rzd with N due to the shell closure at iV = 50 is observed). Therefore, it seems necessary 



to investigate an isospin dependent correction for the scattering of rzd- in Ref- [0,|n| , the 
isospin effect has been considered based on the A 1 / 3 law. However, considering the fact that 
the Z 1 ' 3 dependence can describe the nuclear charge radii much better than the A 1 ' 3 law, 
we take the Z 1 ^ 3 dependence as a more reasonable starting point for describing the isospin 
dependence of nuclear charge radii. 

To confirm that isospin dependent Z 1 ! 3 formula to be developed for nuclear charge radii 
mentioned above is also valid for nuclei far from the /^-stability line, the charge radii in 
exotic nuclei is needed. However, as no such data is available, what we can do is to require 
that our new isospin dependent Z 1 / 3 formula should assort with a reliable and microscopic 
nuclear model. 

The fully self-consistent and microscopic relativistic continuum Hartree-Bogoliubov 



(RCHB) theory, which is an extension of the relativistic mean field (RMF) fll8|-p0"f and 
the Bogoliubov transformation in the coordinate representation |^T|, is a good candidate 
for the present purpose. The RCHB theory can describe satisfactorily the ground state 
properties for nuclei both near and far from the /3-stability line. A remarkable success of 



the RCHB theory is the self-consistent reproduction of the halo in n Li [^2] and prediction 
of the exotic phenomenon - giant halo |23|| . In combination with the Glauber model, the 
RCHB theory successfully reproduces the interaction cross section in Na isotopes p4j and 



the charge changing cross section of C, N, O, F isotopes (ranging from the /3-stability line to 
the neutron drip line) on the target of 12 C at 930 MeV/u [[25|,[2(J. These successes encourage 



us to apply the RCHB theory for the description of charge radii of nuclei both close to 
and far from the /3-stability line and check its validity for the data available and provide 
information for nuclei far away from the stability line. 

The detailed formalism and numerical techniques of the RCHB theory can be found in 
Ref. ETJ and the references therein. In the present calculations, we follow the procedures 



in Ref. [pl| , P3 , M and solve the RCHB equations in a box with the size R = 20 fm and a 



step size of 0.1 fm. The parameter set NL-SH [27] is used, which aims at describing both 



the stable and exotic nuclei. The density dependent 5-iorce in the pairing channel with 
Pq = 0.152 fm -3 is used and its strength Vq is fixed by the Gogny force as in Ref. pi ]. The 
contribution from continua is restricted within a cut-off energy E cut ~ 120MeV. 

As typical examples, we studied the even-even Ca, Ni, Zr, Sn and Pb isotopes ranging 
from the /3-stability line to neutron drip line. The two neutron separation energies S2 n 
is one of the essential quantities to test a nuclear model. In Fig. 2, the calculated S 2n 
(open symbols) of the even-even Ca, Ni, Zr, Sn and Pb isotopes by the RCHB theory are 
compared with the data available (solid symbols) |2J| , where a satisfactory agreement is seen. 
Particularly the deviation between the calculated binding energies with the data available 
is within 1%. In the present calculation, the neutron drip-line nuclei are predicted at 72 Ca, 
98 Ni, 140 Zr and 176 Sn, respectively. In the ^n versus N curve for each isotopic chain, there 
are some kinks due to the neutron shell or subshell closure. For example, the closed shells 
at N = 20, 28 and subshell at N = 40 correspond to kinks in the S , 2 n versus N curve for 
Ca isotopes at 40 Ca, 48 Ca and 60 Ca, respectively. While the kink at iV = 20 for 40 Ca may 
be also due to the Wigner term for N = Z = 20. However, there are no kinks at 70 Ca and 
176 Sn, which indicate the disappearance of magic number N = 50 and 126 for these nuclei 
in RCHB. 

The rms charge radii (r 2 ) 1 / 2 obtained from the RCHB theory (open symbols) and the 
data available (solid symbols) for the even-even Ca, Ni, Zr, Sn and Pb isotopes are given in 
Fig. 3. As it could be seen, the RCHB calculations reproduce the data very well ( within 
1.5%). For a given isotopic chain, an approximate linear N dependence of the calculated 
rms charge radii (r 2 ) 1//2 is clearly seen in Fig. 3, which shows that the variation of (r 2 ) 1 / 2 
for a given isotopic chain deviates from both the the simple Z 1//3 dependence and the simple 
A 1 / 3 law (denoted by dashed lines in Fig. 3). Therefore, a strong isospin dependence of 
nuclear charge radii is necessary for nuclei with extreme N/Z ratio. 

In Fig. 4, the experimental and RCHB predicted rzd — R c /Z 1 ^ for the proton magic 
isotopes are presented as a function of rj = N/Z. It is clearly seen that the coefficient rzd 
increases linearly with 77 (except some deviations due to deformation or shell effect) and the 



slopes are nearly the same for these isotopic chains. The linear 77 ( or isospin Tz = (N — Z)/2 
) dependence of rzd for an isotopic chain may be understood as the effect of the first order 
perturbation correction of nuclear wave function due to an isospin Tz dependent interaction 
p9| . Based on the analysis of data in the middle and upper panels of Fig. 1 and RCHB 
prediction in Fig. 3 and 4, we propose the following isospin dependent Z 1 ^ 3 formula for 
nuclear charge radii: 

R c = a Z 1 ' 3 [1 + b(r) - rf)] , 77 = N/Z, (4) 

where rf is rj = N/Z for the nuclei along the /3-stability line which can be directly extracted 
from the nuclear mass formula |], a = r zd (l + -§^P 2 ), r* Zd = 1.631(11) fm obtained in upper 
right panel of Fig.l, and b = 0.062(9) obtained from the least square fitting. 

The analysis of the available data using Eq. (f|) with r zd and b thus obtained is displayed 
in the lower right panel of Fig. 1. It is found the data are reproduced better by Eq. (^) than 
by Eq. (Q) (the rms deviation is reduced by about 40%). The same has been done for A 1 / 3 
dependence but with less success ( see the lower left panel of Fig. 1 ) . In Ref . [|i~3l|l7| , a similar 
equation for A 1 ! 3 dependence has been used to describe the isospin dependence of charge 
radii and a better agreement has been achieved by fitting both r^ and b simultaneously. 
But then a simple explanation for ta (saturation property) and 6(isospin effect) in Eq. (|J) 
is missing. It is expected that the modified Z 1 ! 3 formula (Eq. (|4])) will become more useful 
with more and more data obtained for the nuclei far from the /3-stability line. 

In summary, we have systematically investigated the nuclear charge radii with A > 
40. It is clearly seen that the Z 1 ' 3 dependence is superior to the A 1 ' 3 law. A delicate 
scattering of data around R c /Z 1 ^ 3 = 1.631 is infered as owing to the isospin effect and a 
linear dependence of R c /Z 1 ' 3 on N/Z ( or (N — Z)/2 ) is found. This inference are well 
supported by the microscopic RCHB calculation conducted for the proton magic Ca, Ni, 
Zr, Sn and Pb isotopes including the exotic nuclei close to the neutron drip line, which 
reproduce S^n and nuclear charge radii data available well. With the linear dependence 
of the coefficient rzd on N/Z (or (N — Z)/2 ) read from data and RCHB theory, a new 



isospin dependent Z 1 ^ 3 formula for the nuclear charge radii is proposed, which improves the 
description of the data available for nuclei near the /3-stability line and could be very useful 
for new data obtained for nuclei far from the /3-stability line. 
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FIG. 1. The nuclear charge radius data for r& in A 1 ' 3 and rz in Z 1 ' 3 law with and without 
isospin dependence, for the details see the text. 
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FIG. 2. Two-neutron separation energies S2n of even Ca, Ni, Zr, Sn, Pb isotopes as a function 
of N, including the data (solid symbols) from Ref.[28] and the RCHB calculation with a <5-force 
(open symbols). 
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FIG. 3. The rms charge radii versus the neutron number N for even-even Ca, Ni, Zr, Sn, Pb 
isotopes. The RCHB calculation with 5-force is represented by open symbols, while the corre- 
sponding data is denoted by solid symbols. The dashed lines represent the predictions by the A 1 ' 3 
law with rAd = 1.228 fm. 
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FIG. 4. The experimental (solid symbols) and RCHB predicted ( dashed lines ) coefficient 
rzd = Rc/Z l f° r the nuclear charge radii as a function of isospin quantity rj = N/Z in even-even 
Ca, Ni, Zr, Sn and Pb isotopes. An asymptotic behavior is drawn as a solid line. 
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